Zero-temperature Glauber dynamics on small-world networks 
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The zero-temperature Glauber dynamics of the ferromagnetic Ising model on small-world net- 
works, rewired from a two-dimensional square lattice, has been studied by numerical simulations. 
For increasing disorder in finite networks, the nonequilibrium dynamics becomes faster, so that the 
ground state is found more likely. For any finite value of the rewiring probability p, the likelihood of 
reaching the ground state goes to zero in the thermodynamic limit, similarly to random networks. 
The spin correlation f(r) is found to decrease with distance as £(r) ~ exp(— r/A), A being a corre- 
lation length scaling with p as A ~ p~ 0,73 . These results are compared with those obtained earlier 
for addition- type small world networks. 

PACS numbers: 64.60.De, 05.50.+q, 05.70.Ln, 89.75.Hc 



I. INTRODUCTION 

In recent years, researchers have accumulated evidence 
that various kinds of complex systems can be described 
in terms of networks or graphs, where nodes play the role 
of system units and edges represent interactions between 
connected pairs of units. Thus, complex networks have 
been used to model several types of real-life systems,and 
to study various processes taking place on them [l|, 0, H, 
13, [5[. In this context, some models of networks have 
been designed to explain empirical data in several fields 
(sociology, economy, biology, technology), as is the case of 
the so-called small-world networks, introduced by Watts 
and Strogatz in 1998 @|. 

This kind of small-world networks are well suited to 
describe properties of systems with underlying topolog- 
ical structure ranging from regular lattices to random 
graphs @, H[ , by changing a single parameter [§] . They 
consist of a regular lattice, in which a fraction p of the 
bonds between nearest-neighbor nodes are replaced by 
new random links, thus creating long-range "shortcuts" 
0, Q . This procedure generates networks in which one 
finds at the same time a local neighborhood (as in regular 
lattices) and some global properties of random graphs, 
as a small mean distance between pairs of sites. This 
short global length scale has been found to be relevant 
for several statistical physical problems on small-world 
networks, such as spread of infections [HI EH, signal 

and information spreading 




Cooperative phenomena in this kind of networks dis- 
play unusual characteristics, associated to their pecu- 
liar topology 0, H3, [H, M, [H HE p. Thus, a 
paramagnetic-ferromagnetic phase transition of mean- 
field type at finite temperature was found for the 
Ising model on small-world networks derived from one- 
dimensional (Id ) Ipll. [pR [29l ]. as well as from 2d and 3d 
regular lattices [30L l3ll] . 

While the single-spin-flip Glauber dynamics at finite 
temperatures reaches the state given by thermodynamic 
equilibrium, the situation at zero temperature is not so 



clear. In fact, it has been shown that nontrivial phe- 
nomena appear, even for regular lattices j32l. l33l. 133. l35l| . 
In d = 1 one reaches the ordered ground state, but for 
higher dimensions the system may get stuck in a frozen 
state with domains of opposite magnetization 35]. The 
Glauber dynamics has been also investigated on vari- 
ous types of complex networks, in particular in small- 
world networks [3(1 H3], Erdos-Renyi random graphs 
[H, SI, HI; and scale-free networks [M HJ SH, E| . In 
these cases, the system may become trapped in a set of 
ordered domains, without reaching a fully ordered state, 
even at finite system size. Analytical calculations, in par- 
ticular, have been carried out to study several character- 
istics of this zero-temperature dynamics on reg ular lat- 
tices [H, |H, |H, Ell and complex networks [ll 11 El • 
A comparison with numerical simulations has shown that 
mean-field-type descriptions may be inappropriate to de- 
scribe this problem [42|. 

Here we consider the zero-temperature Glauber dy- 
namics in small- world networks rewired from a 2d square 
lattice. This allows us to study the evolution of the sys- 
tem from a regular lattice to a random network by chang- 
ing the rewiring probability p. Earlier works considered 
networks in which links were added to the regular lat- 
tice, so that a random network was not approached in 
the limit of large disorder [36|, [37| • We will look for sim- 
ilarities and differences between the spin dynamics on 
these "addition-type" networks and those studied here. 
A question of particular interest is whether for small dis- 
order (p « 1) the system behaves like a random network, 
or on the contrary it keeps some characteristics of the 
regular lattice. 

The paper is organized as follows. In Sec. II we de- 
scribe the model employed in the paper. In Sec. Ill we 
discuss the ordering process, presenting results for the 
fraction of ordered and active runs, as well as the order- 
ing time. In Sec. IV we give results for the active links, 
in Sec. V we present the spin correlation, and in Sec. VI 
we discuss the effect of the initial magnetization on the 
system evolution. The paper closes with a summary in 
Sec. VII. 
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II. MODEL 

To generate our networks, we start from a two- 
dimensional square lattice. Contrary to the Id chain, it is 
known that in 2d lattices the zero-temperature Glauber 
dynamics may get stuck in a frozen state without reach- 
ing the ground state (a ferromagnetic state, with all spins 
parallel), which makes the dynamics nontrivial [34l. l35| . 
Small-world networks were built up according to the 
model of Watts and Strogatz [f| Q, i.e., we considered 
in turn each of the bonds in the starting 2d lattice and 
replaced it with a given probability p by a new connec- 
tion. In this rewiring process, one end of the selected link 
is changed to a new node chosen at random in the whole 
network. We impose the conditions: (i) no two nodes 
can have more than one bond connecting them, (ii) no 
node can be connected by a link to itself, and (iii) each 
node has at least two connections. With this method we 
obtained networks where more than 99.9% of the sites 
were connected in a single component. Moreover, this 
rewiring procedure keeps constant the total number of 
links in the rewired networks, so that we have an average 
degree (k) = 4 irrespective of the rewiring probability 
p. This allows us to study the effect of disorder upon 
the properties of the model, without changing the mean 
connectivity. For networks generated in the present way 
there is a p-dependent crossover size that separates the 
large- and small-world regimes, and the small-world be- 
havior appears for any finite value of p_[0 < p < 1) as 
soon as the network is large enough (45l.l46l| . All networks 
considered here are in the small- world regime [H, [3(| ■ 

We note that other ways of generating small- world net- 
works from regular lattices have been proposed [13, [H[ . 
In particular, instead of rewiring each bond with prob- 
ability p, one can add shortcuts between pairs of sites 
taken at random, without removing bonds from the regu- 
lar lattice. This method turns out to be more convenient 
for analytical calculations, but does not keep constant 
the mean degree (k), which in this case increases with 
p. The zero-temperature Glauber dynamics of the Ising 
model has been studied earlier on these "addition-type" 
small-world networks. Thus, Boyer and Miramontes [361 ] 
studied interface motion and pinning in the limit p< 1, 
whereas Das and Sen [37| considered densely connected 
small-world networks, generated from one-dimensional 
rings. 

From the 2d square lattice, we generated small-world 
networks of different sizes N = L x L. For convenience, 
we will employ later L instead of N to indicate the sys- 
tem size in some plots. The largest networks used here 
included N — 400 x 400 nodes. We considered periodic 
boundary conditions for the starting regular lattice, and 
then the rewiring process was carried out in a way similar 
to that described in Ref. [f|. Each network is character- 
ized by its adjacency matrix Aij, where Aij = 1 if nodes 
i and j are connected by a link, and — otherwise. 
For a given network, we consider an Ising model with spin 
variables Si = ±1 (i = 1, ■■■,N) located on the nodes of 



the network, i.e., we have a Hamiltonian: 

H = -J2AijSiS j . (1) 

i<j 

This means that each edge in the network represents a 
ferromagnetic interaction between spins on the two linked 
nodes. The spin configuration evolves in time, and at a 
given time t, each site i experiences a local field hi(t) due 
to the spins located at its nearest- neighbor nodes: 

hi(t)=J2 A ijSj(t)- (2) 

3 

The zero-temperature Glauber dynamics is then defined 
as follows. At time step t + 1 one selects at random 
a node from a uniform probability distribution, i.e. all 
nodes are equally likely to be selected, irrespective of 
the particular properties of each one, such as its degree 
k. Note that other ways of choosing the nodes, as a 
deterministic sweep, or according to a degree-dependent 
probability distribution, introduce an undesired bias in 
the spin dynamics. Once a node is selected, the value of 
its spin is updated according to the local field, namely: 

C +1 if hi(t) > 
Si(t + 1) = < -1 i£hi(t) < (3) 
I ±1 with probability | if hi(t) = 0. 

We will measure time as the number of attempted up- 
dates per node, N s = t/N (simulation sweeps), so that 
on average each node tries to change its state once per 
simulation sweep. Several variables characterizing the 
considered model have been calculated and averaged for 
different values of p. In general, we have considered 1000 
simulation runs for each rewiring probability p, but in 
some cases we carried out up to 5000 runs to improve 
the precision of our results, in particular to determine the 
average number of active links (see below). We started 
the simulations by assigning to each node a spin +1 or 
— 1 with 50% probability. In Sect. V, however, we con- 
sider the effect of having a starting finite magnetization 
on the evolution of the system. For each run considered, 
the simulation proceeds while the energy of the system 
is changing, and stops when it is constant for a long 
simulation interval. This interval is taken to be equal 
to 2000./V attempted spin updates, i.e., 2000 simulation 
sweeps. This is what we will call in the following the 
infinite-time limit, to distinguish it from the simulation 
times at which the system energy is still changing. 

III. ORDERING PROCESS 

For the zero-temperature Glauber dynamics on the 2d 
square lattice, it is known that the system reaches ei- 
ther a frozen stripe state with probability w 1/3 or the 
ordered ground state with probability « 2/3 [35j. For 
small-world networks we find that the fraction of or- 
dered configurations depends markedly on the rewiring 



3 




1000 



Simulation time (N s ) 



FIG. 1: Fraction of active spin configurations vs. number 
of simulation sweeps, N s , for several values of the rewiring 
probability p. From top to bottom: p = 0.1, 0.3, 0.5, 0.7, and 
1. Data correspond to simulations on networks with 6400 
nodes. 
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FIG. 2: Fraction of simulation runs reaching the ground state 
vs. network size N in the limit of infinite simulation time. 
Each kind of symbols indicates a rewiring probability p. From 
top to bottom: p — 1, 0.7, 0.5, 0.3, and 0.1. Lines are guides 
to the eye. Error bars are on the order of the symbol size. 



probability. We will call active runs those that have not 
reached the ground state at a given simulation time. For 
a given parameter set (p, N) we will denote the fraction 
of active runs by x a . In Fig. 1 we display function 
of the number of simulation sweeps, for several values of 
p. These results were obtained by averaging in each case 
over 1000 spin configurations on networks including 6400 
nodes. From these results we observe that the system 
gets ordered more frequently as the rewiring probability 
increases. In fact, for p — 0.1 we find at large simulation 
times a fraction of active spin configurations larger than 
0.99, i.e. a fraction of ordered configurations less than 
0.01. By comparing with the known result for the reg- 
ular lattice (p = 0), this means that a small fraction of 
rewired links is enough to keep the system with a large 
probability in a disordered configuration. This is in line 
with the result found in Ref. for addition-type small- 
world networks, in the sense that nodes with long-range 
connections strongly affect the motion of interfaces and 
thus inhibit reaching the ground state. 

For larger p the underlying regular lattice is progres- 
sively destroyed, approaching a random network in the 
limit p — > 1. For random networks it is known that the 
long-time behavior of the model is very different from 
that of the 2d square lattice. For such networks the prob- 
ability of the system reaching the ground state goes to 
zero in the large-network limit (N — ► oo) 38]. In Fig. 1 
one also observes that the typical ordering time for small- 
world networks decreases as the rewiring probability is 
raised. This is general for different network sizes, as will 
be discussed below. 

For finite networks, it is interesting to check the 



crossover from the properties of the regular lattice 
to those characterizing the small- world regime. This 
crossover will happen at a rewiring probability p c that 
depends on the system size, and decreases as the system 
size rises (for N — > oo, p c — > 0). Thus, for the system size 
TV = 6400 considered in Fig. 1, one has a change in the 
fraction of active networks at large time, x a , from « 1/3 
for the square lattice (p = 0) to 0.99 for p = 0.1. This 
means that there is a fast change of x a in this region, due 
to the onset of the small-world regime. We have checked 
that, in fact, the crossover appears for this system size 
at p c w 0.02. 

In Fig. 2 we show the fraction of ordered configu- 
rations, 1 function of the network size N, 
in the limit of large simulation time. We display in 
a semilogarithmic plot results for various values of the 
rewiring probability, between p = 0.1 and p — 1. In all 
cases the probability of a spin configuration reaching the 
ground state decreases exponentially as the network size 
increases, and the system remains trapped in spin con- 
figurations with part of the nodes with Si = 1 and the 
rest with Si = — 1. This behavior is similar to that ob- 
served for the Glauber dynamics in random networks at 
T = 0. In fact, Svenson [49( has noticed a freezing in a 
disordered state for Glauber dynamics on random graphs. 
This problem was considered analytically by Haggstrom, 
who showed that the dynamics fails to reach the ordered 
ground state in random networks in the large-size limit 
N — > oo [38| . It is important to emphasize the change in 
x a from the regular lattice to small-world networks in the 
thermodynamic limit. In fact, x a changes from rs 1/3 for 
the square lattice (p = 0) to 1 for any rewiring probability 



4 



400 




50 100 150 

Network size (L) 



FIG. 3: Mean ordering time as a function of network size for 
several values of the rewiring probability p. Symbols indicate 
results derived from zero-temperature simulations. From left 
to right: p = , 0.1, 0.2, 0.3, 0.5, and 1. 
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FIG. 4: Fraction of active links in active runs as a function of 
the number of simulation sweeps. Lines correspond to various 
values of the rewiring probability p. From top to bottom: 
p — 1, 0.7, 0.6, 0.5, 0.1, and 0. These data correspond to 
simulations on networks with 6400 nodes. 



IV. ACTIVE LINKS 



p > 0. The small-world topology thus precludes the sys- 
tem from reaching the ordered ferromagnetic state in the 
large-size limit. The similarity of our results for rewired 
small-world networks with those found for random net- 
works is a difference with addition- type networks, since 
in the latter the connections in the underlying regular 
lattice are preserved for any value of p, and one does not 
recover the behavior of random networks. 

Another interesting quantity is the mean ordering time 
of the system, N . It is defined as the average simulation 
time employed by the system to reach a ground-state con- 
figuration, calculated for the runs that actually become 
ordered. The mean ordering time is displayed in Fig. 3 
as a function of the network size for several values of p, 
between and 1. For a given rewiring probability, N D 
increases linearly with L, i.e. N a ~ y/N. For a given 
network size, the system gets ordered faster for larger p. 
This means that for small p there appear less ordered 
configurations (see Fig. 2), and those that appear take 
longer times to reach the ground state. This is an inter- 
esting difference between the behavior of x a and N a for 
small p. For a given network size, N a evolves in a con- 
tinuous manner from p = to p — 1, as shown in Fig. 3. 
However, the fraction of active runs x a for p = behaves 
in a way different from that for p > 0, since in the former 
case x„ converges to sa 1/3 as the network size increases 
[35l |5Q| , and in the latter x a — > 1 in the large-size limit 
(see above). 



The failure of the system to reach an ordered state 
for system size N — * 00 could be due to the presence of 
small clusters (communities) of nodes tightly connected 
with each other, but loosely linked to the rest of the net- 
work. Such clusters could become ordered independently 
of the rest of the system, thus giving a disordered state, 
and hindering the reach of the ground state. This is the 
mechanism discussed in Ref. [39 in connection with the 
Glauber dynamics on random networks. To address the 
validity of this assumption, it is useful to consider the ac- 
tive links in a network, defined as those connecting nodes 
with opposite values of Si. In the following we will call 
xi the fraction of active links, averaged over active runs. 
The evolution of xi along the simulation runs is presented 
in Fig. 4 for various values of the rewiring probability and 
networks including 6400 nodes. We note first that, for a 
given simulation time N s , the fraction of active links in- 
creases as p is raised. For p < 0.5, x\ decreases as time 
proceeds and eventually converges to a finite value. For 
p > 0.5 we find a decrease in xi at short times, reaching 
a minimum for N s ~ 100, and rising latter to saturate to 
a finite value at large times. In principle, for the order- 
ing process advancing as time proceeds, the number of 
actives links in active configurations decreases, as found 
for p < 0.5. However, a fast decrease in the fraction of 
active runs can cause a rise in xi , since the latter is calcu- 
lated over runs still active, and those getting ordered are 
expected to have few active links. In fact, the minimum 
in xi observed for p > 0.5 occurs at the same simulation 
times as the fast decrease in x a shown in Fig. 1. This 
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FIG. 5: Fraction of active links in active runs vs rewiring 
probability p, in the limit of infinite simulation time. Symbols 
represent results of simulations for networks with different 
sizes: squares, L — 100; circles, L — 300; triangles, L = 400. 
Error bars are on the order of the symbol size. Lines are 
guides to the eye. 



minimum in xi is reminiscent of that found by Castel- 
lano et al. [391 ] for random networks with average degree 
(k) = 7. 

The long-time value of xi is presented in Fig. 5 as a 
function of the rewiring probability p for three different 
system sizes. Symbols indicate results of our simulations: 
squares for L = 100, circles for L = 300, and triangles 
for L = 400. For L = 100 we find that xi increases 
with p, but has a plateau at intermediate values of the 
rewiring probability. This plateau, however, seems to 
be a finite-size effect, since it tends to disappear as the 
network size is increased. In fact, for L — 400, it is 
almost inappreciable; one has a fast rise in the fraction 
of active links xi between p = and p ~ 0.3, and the 
rise is much slower for p > 0.4. Finally, in the large- 
disorder limit (p = 1) xi converges to a value of about 
0.17, irrespective of the system size L. Interestingly, a 
clear finite-size effect appears for intermediate values of 
p, but is very small in both limits p — > (regular lattice) 
and p — > 1 (random networks). Close to p = we find a 
linear increase in xi, with a slope dxi/dp — 1.02 ± 0.02. 

We now go back to the question posed above on the 
presence of small clusters of nodes that could get ordered 
independently of the rest of the network. If this was the 
reason for avoiding the full ordered state, one would ex- 
pect to find a small fraction of active links, which should 
decrease for rising network size. This is not the case, as 
shown by our results displayed in Fig. 5. In line with 
this, we have also observed in our simulations that the 
long-time magnetization goes to zero when increasing the 
system size, as expected for the presence of different large 
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FIG. 6: Spin correlation £(r) vs distance r for active runs on 
small- world networks with 9 x 10 4 nodes, in the limit of large 
simulation time N s — * oo. Error bars are in the order of the 
symbol size. 



ordered regions in the network. This conclusion agrees 
with that given in Ref. [3^ for random networks, but here 
it extends to small-world networks in the whole range 
from p = to 1 . 

In connection with this, it is important the question 
whether the disordered (metastable) state reached after 
long simulation time is a frozen state or not. In our sim- 
ulations we find that the disordered state is active, in 
the sense that some spins continue flipping forever, with- 
out changing the energy of the system. Both spin states 
(+1 and —1) coexist, with several connected domains of 
identical spins. The continuous flipping happens at the 
borders of these domains. This behavior is similar to 
that observed for regular lattices for d > 2 and for ran- 
dom networks [3£|. This means that a small fraction of 
rewired links (p <S 1) is enough to recover the behavior 
of random networks, instead of the freezing observed for 
the Glauber dynamics on the 2d square lattice [35|. 



V. SPIN CORRELATION 

From the results presented above, it is clear that the 
fraction of ordered configurations on small-world net- 
works goes to zero in the large-size limit, and thus the 
ground state is not reached in the thermodynamic limit. 
A direct way of displaying the lack of long-range ferro- 
magnetic ordering is by looking at the spin correlation vs 
distance for various values of the rewiring probability p. 
We define £ as 

£(r) = (SiSJr , (4) 
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FIG. 7: Correlation length A vs rewiring probability in a 
logarithmic plot. Data points were obtained from the decay of 
the spin correlation £(r) with distance, for networks including 
9 x 10 4 nodes, in the limit of large simulation time N s — > oo. 
Error bars are on the order of the symbol size. 
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FIG. 8: Fraction of active links vs network size for several 
values of the initial magnetization. Symbols indicate results 
derived from simulations for a rewiring probability p — 0.1, 
in the limit of infinite simulation time. From top to bottom: 
/+ = 0.5, 0.52, 0.53, 0.55, and 0.6. For /+ > 0.6, x l decreases 
very fast for increasing network size (not shown). 



where the subscript r indicates that the average is taken 
for the ensemble of pairs (i, j) of sites at distance r. Note 
that r = d/d,Q refers here to the dimensionless distance 
between sites in the starting regular lattice, not to the 
actual topological distance or minimum number of links 
between nodes in the rewired networks (do is the distance 
between nearest neighbors). The correlation £(r) in the 
long-time limit is shown in Fig. 6 for several values of the 
rewiring probability p. After a short transient for small 
r, £(r) displays a logarithmic decrease for increasing dis- 
tance. At distances longer than those presented in Fig. 6, 
£(r) saturates to a constant value, which is a finite-size 
effect and should disappear in the limit N — > oo. Note 
that the results displayed in Fig. 6 correspond to val- 
ues of the rewiring probability p <C 1 ■ For larger p, the 
decrease in £(r) as a function of distance is very fast. 

According to the results shown in Fig. 6, we find 
a region where the spin correlation scales as £(r) ~ 
exp(— r/A), with a correlation length A dependent on the 
rewiring probability p. This dependence of A on p is 
shown in Fig. 7 in a logarithmic plot, where one observes 
that the correlation length follows over two decades a 
power law A ~ p~ a , with an exponent a = 0.73 ± 0.02. 

Boyer and Miramontes [36| interpreted the presence 
of different spin domains on their addition-type small- 
world networks in terms of the "influent" nodes. These 
are nodes that have long-range connections, and strongly 
affect the motion of domain interfaces. Moreover, they 
argued that the exponent a should be 2/3, on the ba- 
sis of the random distribution of influent nodes over the 
square lattice. In fact, they found a value a = 0.64 from 
their numerical simulations, close to 2/3. The exponent 



yielded by our simulations on rewired networks is some- 
what higher. This can be understood by taking into ac- 
count that in our case the influent nodes play a role simi- 
lar to that played in addition-type small-world networks, 
but also connections in the underlying lattice are pro- 
gressively destroyed as the rewiring probability increases. 
This affects the correlation between sites in the lattice, 
which is effectively reduced and the exponent a increases 
with respect to that expected for the full square lattice 
(A decreases faster for rising p). 



VI. INFLUENCE OF INITIAL CONDITIONS 

The results shown above were obtained for initial spin 
configurations with zero average magnetization (spins 
with 50% probability for +1 and —1). One can also con- 
sider starting spin configurations with different probabil- 
ities for spins +1 and — 1, and study the evolution of the 
system under the Glauber dynamics at T = 0. A study of 
the effect of initial conditions on this type of dynamics in 
various kinds of complex networks has been carried out 
earlier by Uchida and Shirayama [40j . These authors em- 
phasized that both the initial conditions and the network 
structure are equally relevant to determine the evolution 
of the system. This evolution may be nontrivial, mainly 
in the presence of nodes with large degree, as happens in 
scale-free networks. 

We now consider initial configurations on our small- 
world networks, where the probability /+ of Si = +1 is 
higher than 0.5, and study their evolution with time. In 
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FIG. 9: Average magnetization M of spin configurations in 
active runs as a function of network size, for various values of 
the initial magnetization. Symbols show results derived from 
simulations on networks with a rewiring probability p — 0.1, 
in the limit of infinite simulation time. From top to bottom: 
/+ = 0.6, 0.55, 0.53, 0.52, and 0.5. For /+ > 0.6, M converges 
fast to 1 for increasing network size (not shown). Dashed lines 
are guides to the eye. 



Fig. 8 we present the fraction of active links as a function 
of the network size for several values of / + between 0.5 
and 0.6. These results were obtained for a rewiring prob- 
ability p = 0.1 and in the long-time limit N s — > oo. For 
/ + = 0.5 (random initial configuration), the fraction Xi 
converges to a finite value as the network size increases, 
as shown above. For / + > 0.5, however, xi decreases 
exponentially as N is raised, indicating that the system 
evolves to a single ordered domain, eventually converging 
to the ground state. 

To check this point we have calculated the average 
magnetization M over active networks as a function of 
the system size. We define M as M = (|M|) ac , with 
M = J^i ln Fig- 9 we show M for the same rewiring 
probability and /+ values presented in Fig. 8. Again the 
results correspond to the long-time limit. For random ini- 
tial conditions the average magnetization decreases and 
goes to zero for N — > oo, as indicated above (for L = 300, 
we find M = 0.057, not shown in the figure). On the con- 
trary, for /+ > 0.5, the average magnetization increases 
as the network size is raised, converging to the value ex- 
pected for the fully ordered system (M = 1). This is 
consistent with the decrease observed in the fraction of 
active bonds xi for rising N, shown in Fig. 8. More- 
over, this behavior is reminiscent of that observed for 
random networks in Ref. in the sense that a small 
finite value of the starting magnetization is enough to 
drive the system to an ordered state. Such an evolution 
can be different for other kinds of complex networks • 

Uchida and Shirayama 40] found that, depending on 



the type of network, the final dynamics may reach a 
metastable state involving two coexistent spin states 
with several connected domains of identical spins, the 
marginal vertices of which flip continuously. This is in 
fact the behavior found for our small-world networks and 
presented above for random initial conditions (/+ = 0.5) 
in section IV. For / + > 0.5, however, we have found that 
the system converges to an ordered state in the thermo- 
dynamic limit N — > oo [M — > 1), but has a finite proba- 
bility of being disordered for any finite size. Both Figs. 8 
and 9 indicate that, for a given value /+ > 0.5, the ap- 
pearance of a single spin domain is favored as the system 
size increases. However, the convergence to the ordered 
state for increasing system size is fast for p > 0.6, but 
becomes very slow as one approaches the random initial 
conditions, and finally the system remains in a disordered 
state for / + = 0.5 and N — > oo (see above). Note that for 
f + > 0.5 the behavior of the system is opposite to that 
shown above for /+ = 0.5, in the sense that for random 
initial conditions the system reaches the ground state less 
frequently for increasing size N (see also the increase in 
Xi shown in Fig. 8 for /+ = 0.5 and decrease in M in 
Fig. 9, as rises). 



VII. SUMMARY 

We have studied numerically the Glauber dynamics of 
the ferromagnetic Ising model on small- world networks, 
generated by rewiring links in a two-dimensional square 
lattice. We have found that the behavior of the model 
departs from that known for the regular lattice, even for 
a small fraction of rewired links. In fact, for any rewiring 
probability p, the fraction of ordered runs for large time 
(N s — > oo) and large system size (N — > oo) goes to zero, 
contrary to the regular lattice, where this fraction con- 
verges to a finite value (« 2/3). For finite networks, the 
system gets ordered more frequently as the rewiring prob- 
ability increases and one approaches a random network. 

The spin correlation on the underlying lattice is found 
to decrease as £(r) ~ exp(— r/A), with a correlation 
length that depends on the rewiring probability as A ~ 
p~°' 73 . The exponent giving the dependence of A on 
p is close to but different from that derived earlier for 
addition-type small-world networks. 

We have analyzed the influence of the initial condi- 
tions on the behavior of the system. We found that for 
any small deviation from random initial conditions, the 
system evolves to the ordered state in the limit of large 
system size and simulation time, contrary to the behavior 
obtained from random initial conditions. 

The apparently simple Glauber dynamics at T = 
displays a nontrivial behavior, that depends upon the 
type of network on which it is defined. This is true even 
for regular lattices with d > 1. For small- world net- 
works rewired from 2d lattices, we find a behavior of the 
Glauber dynamics similar to that corresponding to ran- 
dom networks. This refers in particular to the nature 
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of the disordered state obtained at long times, which is 
characterized by the presence of ordered regions, with 
spins at the interfaces between domains flipping forever 
without changing the energy. This behavior is differ- 
ent from the frozen state reached for 2d regular lattices, 
but is reminiscent of that for lattices with dimension- 
ality larger than 2. An important result is that for our 
small- world networks, a small fraction of disorder (in fact 
any rewiring probability p > 0) is enough to recover the 
behavior known for random networks. 

For small-world networks generated from lattices with 
a dimensionality different from 2, we expect a behavior 
qualitatively similar to that presented here for d = 2. 
However, the actual details of the nonequilibrium order- 
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